Self— similar and charged spheres 
in the free— streaming approximation 

W. BarretcQ and L. Rosale^l 



Abstract 

We evolve nonadiabatic charged spherical distributions of matter. Dis- 
sipation is described by the free-streaming approximation. We match a 
self-similar interior solution with the Reissner-Nordstrom-Vaidya exte- 
rior solution. The transport mechanism is decisive to the fate of the 
gravitational collapse. Almost a half of the total initial mass is radiated 
away. The transport mechanism determines the way in which the electric 
charge is redistributed. 
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1 Introduction 

There is a renewed interest in the study of self-gravitating spherically symmet- 
ric charged fluid distributions (pQ, [2] and references therein). In self-gravitating 
systems the electric charge is believed to be constrained by the fact that the 
resulting electric field should not exceed the critical field for pair creation, 
10 16 Vcm _1 [3J. This restriction in the critical field has been questioned [1]- 
[7] and does not apply to phases of intense dynamical activity with time scales 
of the order of (or even smaller than) the hydrostatic time scale, and for which 
the quasistatic approximation [8] is clearly not reliable as in the collapse of very 
massive stars or the quick collapse phase preceding neutron star formation. 
Electric charge has been studied mostly under static conditions [SJ— [TT] . Of 
recent interest are charged quasi-black holes [12] and the electrically charged 
extension to quasispherical realization [T3J. Distributions electrically charged 
can be considered in practice as anisotropic [14], [15]. Some authors combine 
anisotropy and electric charge p], [16], [17] but not as a single entity by means 
of an equation of state. 

The electric field has been postulated to be very high in strange stars with 
quark matter [TH] , [T5] , although other authors suggest that strange stars would 
not need a large electrical field [20] . The effects of dissipation, in both limiting 
cases of radiative transport, within the context of the quasistatic approximation, 
have been studied in [3T] . Using this approximation is very sensible because the 
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hydrostatic time scale is very small, compared with stellar lifetimes, for many 
phases of the life of a star. It is of the order of 27 minutes for the sun, 4.5 
seconds for a white dwarf, and 104 seconds for a neutron star of one solar mass 
and 10 km radius [22], [25]. However, such an approximation does not apply to 
the very dynamic phases mentioned before. In those cases it is mandatory to 
take into account terms which describe departure from equilibrium, i.e., a full 
dynamic description has to be used [24"j . 

We are concerned in this paper with configurations out of (static) equilibrium 
with intense dynamical activity. We use the approach of considering the electri- 
cally charged matter distribution as an anisotropic fluid. It is well known that 
different energy-momentum tensors can lead to the same spacetime [25]-|29]. 
Under spherical symmetry, electric charge can be considered as a special case of 
anisotropy |14j , [2] . To obtain a dynamical model we consider the free-streaming 
approximation as the transport mechanism, and a self-similar spacetime for the 
inner region. We explore the fate of the gravitational collapse. Since dissipa- 
tion is introduced by the free-streaming of radiation and local anisotropy by 
the electric charge, it is not necessary to consider in our context any hyperbolic 
theory of dissipation as for heat flow or viscosity. In many circumstances the 
mean free path of particles transporting energy may be large enough to justify 
the free-streaming approximation. We do a comparison with previous work [14] 
in which the transport mechanism was diffusive. 

Is it well known that the Einstein field equations admit nomothetic motion 
[3"0]-|32j. Applications of self-similarity range from modeling black holes to 
producing counterexamples to the cosmic censorship conjecture [33]-[44]. It is 
well established that in the critical gravitational collapse of an scalar field the 
spacetime can be self-similar [I5]-[I7]. We have applied characteristic methods 
to study the self-similar collapse of spherical matter and charged distributions 
05]— [50], [2]. The assumption of self-similarity reduces the problem to a sys- 
tem of ODE's, subject to boundary conditions determined by matching to an 
exterior Reissner-Nordstrom-Vaidya solution. Heat flow in the internal fluid is 
balanced at the surface by the Vaidya radiation. One simulation [50] illustrates 
how a nonzero total charge can halt gravitational collapse and produce a final 
stable equilibrium. It is interesting that the pressure vanishes in the final equi- 
librium state so that hydrostatic support is completely supplied by Coulomb 
repulsion. Another possible final state is extremely compact and oscillatory 
with non zero pressure 14]. In this last case electric charge redistribution in 
the fluid is possible. We explore here if these results depend on the mechanism 
of transport. 

This work is organized as follows. In Sec. 2 we write the field equations for 
electrically charged interior fluids as seen by Bondian observers. We also show 
in this section the junction conditions with the exterior spacetime and sketch 
the general procedure to get physical variables. In Sec. 3 the set of surface 
equations are presented for parametrized self-similar solutions. Modeling is 
performed in section 4 to discuss results and conclude in section 5. 
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2 Field equations 



We proceed now to describe the matter distribution, the inner and outer line 
elements and the field equations. 



2.1 Interior spacetime 



Our starting point is the Bondi approach to study the evolution of gravitat- 
ing spheres |51j . Let us consider a nonstatic distribution of matter which is 
spherically symmetric. In radiation coordinates [52] the metric takes the form 



e 2fi [{V/r)du 2 + 2dudr] - r 2 (d8 2 + sin 2 9d(/) 2 ), 



(1) 



where j3 and V are functions of u and r. Here, 
I + , with u =const. on the outgoing cones, r = a 



u — x° is the Bondi time at 
1 is the surface-area distance 



and 



,3 ; 



is the usual polar coordinates. We use 



geometrized units (c = G = 1). To get physical input we use Bondian observers 
[53"] , introducing the Minkowski coordinates (T, X, Y, Z) by 



dT 
dX 
dY 
dZ 



e f} [(V/r) 1/2 du+ (r/V) 1/2 dr], 

e P(r/V) 1/2 dr, 

rd6. 



(2) 
(3) 
(4) 
(5) 



Next one assumes that for an observer moving relative to these coordinates with 
velocity uj in the radial direction, the space contains a charged fluid with energy 
density p, pressure p, electric energy density jx and unpolarized energy density 
e traveling in the radial direction. For this comoving observer, the covariant 
energy tensor is 



/ p + fi + e -e \ 

-e p-fi + e 

p + fi 

\ p + fi J 



(G) 



where /t = E 2 /8tt and E — Q/r 2 is the local electric field being Q(u,r) inter- 
preted naturally as the charge within the radius r at time u, which satisfies the 
conservation equation 

dv 

DQ = Q.u + —Q.r = 0, (7) 
du 

where the a comma represents partial derivative respect to the indicated coor- 
dinate, and the matter velocity is given by 



dr 
du 



V uj 
r 1 — uj 



(8) 
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Observe that the adiabatic fluid has a diagonal covariant energy tensor (p,p r ,Pt,Pt) 
with p = p + fi, p r = p — jl and pt = p + p\, which is exactly the same as for 
an anisotropic fluid [54]. Clearly the electric charge produces local anisotropy, 
contributing to the matter energy density and pressure. If electric charge is zero 
we recover the Pascalian character of neutral matter, that is, isotropy. 
Now, we define the mass function as 

rh= l -{r~Ve~^), (9) 

related with the usual total mass by means of [3"]. [55] 

Q 2 

m,T = m + — — , (10) 
2r 

which is the generalization of the Misner-Sharp mass for the charged case pQ. 
Thus, the field equations can be written as fPI].[5"6"] 

p+p r u 2 „ (l+a;) _ e~ 2 ^rh iU m r 
1-uj 2 £ (l-w) _ 47rr(r-2m) Anr 2 ' [ ' 

rh T , , 

*=4^' (12) 



Pt= - ^-P,u r e- 20 + -^(l-2m/r)(2/3, rr + 4/3 2 r - /3, r /r) 

47T 07T 

+ [3/3 r (l — 2m r ) — m rr ], 

8irr 



(14) 



where 



and 



Pr - up 

l+LU ' 



LUp r 



1 



(15) 



(16) 



Observe that the field equations ( 11 )-( 14 1 are exactly the same as for anisotropic 



matter [55] , [5T] , that is, electric charge can be interpreted as local anisotropy of 
the fluid. That is possible because of the mass definition given by ( 27 ) , otherwise 
electric charge has to be viewed as part of the metric. From this point of 
view, electric charge is formally an additional physical variable which contributes 
clearly to the matter energy density and pressure. Also this interpretation is 
better understood physically due to the Bondi point of view about the comoving 
reference system. This procedure, except for the mass definition, was the same 
followed to interpret viscosity as anisotropy [56] . 
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2.2 The exterior spacetime and junction conditions 



We consider that the spherically symmetric distribution of collapsing charged 
fluid is bounded by the surface E. Outside E we have the Reissner-Nordstrom- 
Vaidya spacetime, that is, all outgoing radiation is massless, described by 

ds 2 = (1 - 2M(u) /r + q 2 /r 2 )du 2 + 2dudr 

- r 2 (d6 2 + sin 2 6 d<p 2 ), (17) 

where M(u) and q denote the total mass and charge, respectively. 

It can be shown that the junction conditions for the smooth matching of ([I]) 



and (17) on E implies [62] 



p = 0, Q = q, m(u, r) = M(u), /3 = 0, (18) 

where = means that both sides of equation are evaluated on E. It is remark- 
able that the zero pressure on E is equivalent to the continuity of the second 
differential form equation 

+ (i- + £ s* e»> 

Up to this point, within the Bondi framework and spherical symmetry, all 
the written equations are general. We have five physical variables (p, p, u>, e, Q) 
and two metric functions (m, (3), for which we have the five equations ([7]), ( |11| 



(14). Thus, additional assumptions are necessary to solve the characteristic 
initial value problem. In the next section we suppose that the spacetime is self- 
similar to illustrate how the approach works and explore the influence of the 
transport mechanism. But before any assumption to get the metric functions 
we sketch the general procedure to obtain physical variables. 



2.3 Getting physical variables 

The algorithm to calculate the physical variables once we get the dynamics at 
E and the metric functions everywhere is as follows: 

1. Specifying an initial charge distribution (any) we calculate all physical 
variables (p, p, u, i) at any piece of material; 

2. We integrate numerically equation ([T]) to advance Q in time; 

3. Once again we get all physical variables up to the integrate time at any 
piece of material. 
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3 Self— similarity and surface equations 



Self-similarity is invariably denned by the existence of a homothetic Killing 
vector field 30 . A homothetic vector field on the manifold is one that satisfies 
£ £g =2ng on a local chart, where n is a constant on the manifold and £ denotes 
the Lie derivative operator. If n ^ we have a proper homothetic vector field 
and it can always be scaled to have n = 1; if n — then £ is a Killing vector 
on the manifold [571 [58J [5§] . So, for a constant rescaling, £ satisfies £^g =2g 
and has the form £ = A(u, r)d u + X(u, r)d r . If the matter field is a perfect fluid, 
the only equation of state consistent with £^g =2g is a barotropic one [30] . 
The homothetic equations reduce to £(X) — 0, £(Y") = 0, A = r and A = A(u), 
where X = rh/r and Y = Ae 213 /r. Therefore, X = X{Q) and Y = Y{Q are 
solutions if the self-similar variable is defined as £ = r exp(— J du/A). Here we 
assume that X = Ci£ fc and that Y = C^C, where C\, C%, k and I are constants. 
This power-law dependence on £ is based on the fact that any function of £ is 
solution of £^g =2g. Demanding continuity of the first fundamental form we 
get the following metric solutions [50], [T3]: 



fe+i 

r 1 



m = m s ^— J (20) 
and 

/ \ 1+1 

• w 

where k and I are constants; the subscript indicates that the quantity is evalu- 
ated at the surface E. Thus r^(u) represents the radius of the distribution. 
Clearly the metric functions and consequently the matter variables and electric 
charge are determined up to the time dependent functions ms and rs for which 



we have the surface equations obtained from (11) and fl8| evaluated at E 



m E = -47rr|e E (l-^V^ ( 22 ) 



s 



r s J ' ' " 2r 2 

and 



^=(1-^)^-. (23) 
V r E J 1 - 

where = 6^(1 + wa)/(l -ws). Using the dimensionless variables 



m s (0)' ms(0)' 
including the dimensionless time and electric total charge 

U=— ; C=-^ , 

m s (Q)' m E (0)' 

the last two equations can be written as 



dM T ( 2M\ dR C 2 

dU = - L [ 1 -^) + dU2W (24) 
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Figure 1: Evolution of the radius R for k = 0.35, I = 0.5, C = 0.292 and 
R(0) = 2.923. 
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Figure 2: Evolution of the dimensionless energy density p — m s (0) 2 /5 (multiplied 
by 10 2 ) for k = 0.35, I = 0.5, C = 0.292 and R(0) = 2.923. 
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Figure 3: Evolution of the dimcnsionless radial pressure p = ms(0) 2 j5 (multi- 
plied by 10 2 ) for k = 0.35, I = 0.5, C = 0.292 and R(0) = 2.923. 
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Figure 4: Evolution of the dimensionless energy flux e = m s (0) 2 e (multiplied 
by 10 2 ) for k = 0.35, I = 0.5, C* = 0.292 and R(0) = 2.923. 



10 



■0.02 
■0.04 
■0.06 - 
■0.08 - 
■0.10 
■0.12 h 
■0.14 
■0.16 
■0.18 
■0.20 h 
■0.22 
■0.24 



r/r 2 =0 . 6 
r/r 2 =0 . 7 
r/r 2 =0 . 8 
r/r 2 =0.9 
r/r y =l .0 



10 20 30 40 

U 



50 60 70 80 



Figure 5: Evolution of the radial velocity oj for k = 0.35, I = 0.5, C = 0.292 
and R(0) = 2.923. 
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Figure 6: Evolution of the dimensionless charge function Q — m s (0)Q 
k = 0.35, I = 0.5, C = 0.292 and R(0) = 2.923. 
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and 



dR _ ( 2M 
dU ~ \ ~ Jl 



1 - we' 



where L — 47rrg.es- From (19) we obtain 



w E = 1 - 



2(J + 1)(1-2M/.R) 
2(k + l)M/R~C 2 /R 2 



(25) 



(26) 



Now, from the held equation ( 14 1 evaluated at £ we obtain M as a function of 
i?, 



M 



where £ = (/ + l) 2 and ^ = 2£ • 



1 

fc(4 



£i? 

3^- 



C 2 
I? 



(27) 



( 22 ) and rt23 



rc 2 


'i 


1" 




\i? 2 


.2 




tp J 1 - w s 



/c) and consequently we get L from 



(28) 



Thus, we have one independent surface equation to integrate numerically if 



we specify initially in some way M or R. We have preference for ( 23 1 , taking 



into account (26) and (27) 



We do not have, a priori, any restriction on values of k, I and C. Only 
physical criteria and expectations on the foregoing models determine our choices, 
as we illustrate in the following Section. For instance, any choice of values for 
k, I and C have to be consistent with — 1 < u < 1, m > 0, p > 0, p > p for any 
light cone. 



4 Modeling 

As a first step for modeling we want to know the minimal radius for a given set 
of (k, I, C) imposing the conditions on the surface of the distribution. It is easy 
to show that ps > imposes the minimum radius: 



£ (2(fc + l) 

if > — 1. It is also important to note that any model based on self-similar 



solutions (20)-(21) is singular at r = 0. We specify an initial profile for the 



dimensionless electric charge function as 

Q(Q) = c(j-) , (29) 

where V is a free parameter. 

We do a general survey for k, I and C . There are values for which the 
distribution apparently is static. For instance, for k — 0.35, I ~ 0.575 and 
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C = 1(T 3 , the sphere has w s « 1.192 x 1(T 7 and L -4.172 x 1(T 8 , and 
stays there indefinitely. The situation remains the same for other values of total 
electric charge. 

We opt for a general model with 

k = 0.35; I = 0.5; C = 0.292; R(0) = 2.923. 

In this case M(0) f=s 1 and Ws(0) s» —0.037. Thus, we integrate numerically the 



surface equation (23 1 using the fourth order Runge-Kutta method, constrained 
by (26) and (27 1. Choosing the free parameter V = 2 we integrate numerically 
the equation l7| using finite differences. The procedure is straightforward and 
standard (see [2] and references therein). It basically consists of using the Lax 
method (with the appropriate Courant-Friedrichs-Levy (CFL) condition). The 
conservation equation dynamics is restricted by the surface evolution. Once 
the charge function is advanced in time we can get any other physical variable. 
Thus the evolution of the whole distribution proceed up to the final time. For 
other choice of parameters and initial conditions the procedure is the same and 
the results displayed in figures 1-6 are representative. 



5 Discussion 

In this paper we consider electrically charged matter as anisotropic matter. 
We explore a dynamical model under the free streaming approximation and 
self similarity within the source. The example we have shown is representative 
of many others varying the initial condition for the radius R, the self-similar 
parameters k, I, the total charge C. Only regions 0.6 < r/rs < 1 satisfy physical 
conditions pointed out by end of Section 3. 

The main motivation for this work were previous results using the same sys- 
tem and solutions but a different transport mechanism, that is, the diffusion 
limit [14]. Heat flow makes the distribution evolve in a very different way: i) 
the electric charge halts the collapse; ii) the distribution becomes dust asymp- 
totically in one special case; iii) the final state is oscillatory; iv) the electric 
charge is redistributed. 

In the present case the electric charge does not change the fate of the gravi- 
tational collapse. The distribution evolves radiating a huge quantity of its mass 
(« 45%) reaching relativistic velocities of collapse. In fact, inner regions are out 
of the physical domain. Although the evolution looks catastrophic no evidence 
of black hole formation appears during the monitored evolution. 

The only common feature for both transport mechanisms is the electric 
charge redistribution, but in an opposite manner as we explain below. 

We can read from figure 6 that for any time the interior electric charge, for 
any comoving space marker r/r-£, is always less than the total electric charge 
enclosed by the boundary surface. Therefore, the electric charge for inner regions 
can in fact increase (or decrease) by means a redistribution, conserving the total 



14 




Figure 7: Dimensionless charge function Q = ms(0)Q for k — 0.35, I — 0.5, 
C = 0.292 and R(0) = 2.923, as a function of comoving markers r/r s for initial 
and final time U. The electric charge redistributes in the interior space with 
time. 
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electric charge. In general, the electric charge grows from zero to the maximum 
value (the total) in all cases and, in the present (jAs) 2 initially and near 

(r/rs) 4 / 5 later, conserving partially and totally the electric charge, as dictated 
by Eq. (7). It is also important to consider that redistribution occurs while the 
sphere is contracting. Figure 7 shows clearly these issues. In reference [5] an 
analogous behavior is shown in other context, using other method of solution 
but with the same transport mechanism (see figures 15 and 16 there). 

Why the differences found in the redistribution of charge depend on the 
transport mechanism? Clearly, the profiles of the flux energy density for the 
present case (figure 4) and the heat flux as reported in [TJ] (figure 4 there), are 
responsible for opposite radial velocity profiles (figure 5 in [U and figure 5 here) 
and in consequence opposite electric charge redistribution (figure 6 in [Tl] and 
figure 6 here), conserving in all cases the electric charge. We mean by opposite 
that if the free streaming (heat flow) is lesser (greater) at the interior, the radial 
velocity is greater (negatively) at the interior, showing the opposite behavior for 
the diffusion mechanism, that is, lesser (negatively), all this for time-windowing 
where comparison applies. Thus, opposite electric charge redistribution means 
that between fixed extremes the inner electric charge per shell increases with 
evolution for the streaming out and diminishes for the diffusion transport mech- 
anism. 

Generic numerical solutions for homothetic motion are currently under con- 
sideration and will be reported elsewhere. 
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